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Abstract. We present two calculations for a class of robust homoclinic cycles with 
symmetry Z„ K , for which the sufficient conditions for asymptotic stability given 
by Krupa and Melbourne are not optimal. 

Firstly, we compute optimal conditions for asymptotic stability using transition 
matrix techniques which make explicit use of the geometry of the group action. 

Secondly, through an explicit computation of the global parts of the Poincare map 
near the cycle we show that, generically, the resonance bifurcations from the cycles 
are supercritical: a unique branch of asymptotically stable period orbits emerges from 
the resonance bifurcation and exists for coefficient values where the cycle has lost 
stability. This calculation is the first to explicitly compute the criticality of a resonance 
bifurcation, and answers a conjecture of Field and Swift in a particular limiting case. 
Moreover, we are able to obtain an asymptotically-correct analytic expression for the 
period of the bifurcating orbit, with no adjustable parameters, which has not proved 
possible previously. We show that the asymptotic analysis compares very favourably 
with numerical results. 
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1. Introduction 



Heteroclinic orbits between saddle-like invariant sets are of great interest in dynamical 
systems since they generate many kinds of non-trivial behaviour, including intermittency 
and chaotic dynamics fl]. In continuous-time dynamical systems (i.e. ordinary 
differential equations) a non-transversal intersection may give rise to intermittent 
dynamics in the sense that trajectories approach neighbourhoods of the saddle-like 
invariant sets and spent a substantial amount of time there before moving rapidly to a 
neighbourhood of the next invariant set. 

§ Corresponding author (c.postlethwaite@math.auckland.ac.nz) 
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It is well known that, while heteroclinic orbits in generic systems of nonlinear 
ordinary differential equations (ODEs) are of codimension at least one, in sets of ODEs 
containing invariant subspaces they can exist for open sets of parameter values, that is, 
they are codimension zero, and hence are referred to as 'robust' [Ullj. Three situations 
in which such invariant subspaces can arise are (i) due to equivariance with respect to 
a symmetry group [H H], (ii) modelling assumptions such as the permanence of death 
in Lotka-Volterra-type models of population dynamics O [6|, [7], and (iii) structural 
restrictions such as the coupled cell structures investigated recently by Aguiar et al. [8]. 
Detailed results on the existence of homoclinic cycles in particular families of equivariant 
vector fields have been given by many authors, for example P, [ini [H] ■ 

A number of co dimension-one bifurcations have been identified in which robust 
heteroclinic cycles are created or destroyed, or in which their stability changes. Issues 
of stability turn out to be more subtle than might be at first thought. In particular, 
weaker definitions of stability than asymptotic stability turn out to be useful, and indeed 
for some systems describe the generic case. The most prevalent of these weaker notions 
is 'essential asymptotic stability', introduced by Melbourne [12] who gave an example of 
an essentially asymptotically stable robust heteroclinic cycle. Recent work by Driesse 
and Homburg [13] discusses examples of essentially asymptotically stable homoclinic 
cycles which are produced in codimension-one bifurcations from asymptotically stable 
robust homoclinic cycles as a 'transverse' eigenvalue (which we shall define in section [2]) 
crosses the imaginary axis. 

In this paper we discuss a particularly simple class of homoclinic cycles in M". 
Although necessarily restrictive, this class contains the essence of the Lotka-Volterra 
type examples discussed at length by Rabinovich and co-workers jTH [151 [IS [12] as 
models for neural decision-making processes. The class of homoclinic cycles that we 
consider can undergo both transverse and resonance bifurcations and we focus here 
on the resonant case, beginning by proving necessary and sufficient conditions for 
asymptotic stability, improving on the general result given by Krupa and Melbourne [3l 
H]. Such a resonance bifurcation is shown to yield a single periodic orbit which lies 
close to the (now unstable) homoclinic cycle. The stability and direction of bifurcation 
of this periodic orbit depends on coefficients in the Poincare return map which come 
from the global maps, that is, those which describe how trajectories near the homoclinic 
cycle behave outside neighbourhoods of the equilibria on the cycle |I8]. Usually these 
coefficients are impossible to compute analytically, but in this case the calculation 
turns out to be tractable, asymptotically in the limit where the sum of the transverse 
eigenvalues is small. To the best of our knowledge no calculation along these lines has 
been attempted previously. The result of the calculation is to show that the resonance 
bifurcation for this class of simple homoclinic cycles is always supercritical, at least 
when the sum of the transverse eigenvalues is sufficiently small compared to the leading 
expanding and contracting eigenvalues. 

Field and Swift [I9] study in detail a particular example from the class of cycles we 
consider, in the case n = 4. They conjecture that the resonant homoclinic bifurcation 
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is always supercritical. Our results prove this conjecture for an open set of parameter 
values. 

The outline of the paper is as follows. In section [2] we introduce our notation and 
the class of robust homoclinic cycles that we study. We state our result on asymptotic 
stability of such cycles (Theorem 1). Section [3] contains the proof of Theorem 1. In 
section m we present the calculation of the coefficients in the Poincare return map. The 
calculations are reasonably straightforward to follow but become remarkably lengthy. 
Some more detailed parts of the calculations are relegated to the Appendix. In section [5] 
we use our return map calculations to provide an estimate of the period of the bifurcating 
periodic orbit. We also show that the agreement between the analytic leading-order 
estimate and numerical results is very good. We conclude in section [61 

2. Robust homoclinic cycles 

We consider continuous time dynamical systems in the form of F-equivariant ODEs: 

X = /(x), X G M" (1) 

where F C 0{n) is a finite Lie group. An equilibrium ^ G M" of ([1]) satisfies /(^) = 0. 
We begin by giving a number of definitions; these are all standard in the literature, see 
for example [H H] . 

Definition 1 (j)j{t) is a heteroclinic orbit between two equilibria and ^j^i of (J\) if 
(f)j{t) is a solution of (Qp which is backward asymptotic to C,j and forward asymptotic to 

0+1 • 

Definition 2 A heteroclinic cycle is an invariant set X C consisting of the union 
of a set of equilibria {^i, ^m} and orbits {0i, (pm}, where (pj is a heteroclinic orbit 
between and O+i/ '^nd ^m+i =0- require that m > 2. 

In the case m = 1, that is, ^2 = ^i, we say that 0i is a homoclinic orbit to ^1. 

Definition 3 A heteroclinic cycle is a homoclinic cycle if there exists 7 G F such that 
70 = 0+1 for allj- 

Definition 4 For x G MJ^ the isotropy subgroup is 

= {(T G F : ax = x}. (2) 

Definition 5 For S an isotropy subgroup of A, the fixed-point subspace Fix E is 

Fix S = {x G M" : (Tx = X V(T G S}. (3) 

Definition 6 A heteroclinic cycle X is robust if for each j , 1 < j < m, there exists a 
fixed-point subspace, Pj = Fix Sj where C A and 

(i) C,j is a saddle and ^j+i is a sink for the flow restricted to Pj, 
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Table 1. Classification of eigenvalues. P Q L denotes the orthogonal complement in 



P of the subspace 


L. 




Eigenvalue class 


Subspace 


Radial (r) 






Contracting (c) 






Expanding (e) 


V,{e) 


= Pj e Lj 


Transverse (s) 







(a) there is a heteroclinic connection from C,j to ^j^i contained in Pj. 

Importantly, robust heteroclinic cycles may occur as codimension-zero phenomena 
in systems with symmetry. That is, they may exist for open sets of parameter values. 
We define Lj = Pj^i H Pj and clearly ^j G Lj. Following [3], the eigenvalues of the 
linearisation of f{x) about each equilibrium can be classified according to the subspaces 
in which the eigenspaces lie, as shown in table [H 

The specific differential equations we consider in this paper are of the form: 

X = f{x) = + Q{x), xeM" (4) 

where /i G M+ (and so, following a rescaling of time, can be set equal to 1), n > 3, and 
Q{x) is a F-equivariant polynomial, that is, 

lQ{x) = Qicfx), 7 e F, 
which contains only nonlinear terms. The group F is of the form 

F = Z„ K A„, (5) 
where A„ = Z2 acts on M" as n reflections Kj, j = 1, . . . , n: 

Kj , . . . , Xj , . . . , Xn) (3^1 ) . . . ) Xj , . . . , Xfi) , (6) 

The actions of the reflections Kj mean that every coordinate hyperplane is a fixed-point 
subspace, and hence invariant under (jll). The Z„ subgroup of F is generated by a cyclic 
permutation element p which acts as: 

. . . , 2^n) (2^71; Xi, • • • 5 Xn—l)' (7) 

It is clearly sufficient to consider the dynamics restricted to the domain 

M!^ = {(a;i,...,x„) GM"|xi,...,a;n>0}. 

We now describe the robust homoclinic cycles X that we study in this paper. 
For now, we suppose that it is possible to construct a vector field with the following 
properties. In section H] we give an explicit example of ODEs that contain such a cycle. 
We make the following assumptions on the ODEs (jlj): 

(HI) There exist n equilibria ^j, j = 1, . . .n, related by the symmetry element p: 

0+1 = Pl- 
each with isotropy exactly A„_;, (1 < / < n — 2), and therefore with / non-zero 
coordinates. 
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(H2) The unstable manifold of ^j, is one-dimensional, and has isotropy A„_i_x, 

that is, points on the manifold have (Z + 1) non-zero coordinates. 

Without loss of generality, we can choose the non-zero coordinates of to be 
Xi, . . . and the non-zero coordinates of W^{^i) to be Xi, . . . , x^+i. Therefore and 
^2 are contained in the subspace Pi, spanned by We make the further 

assumption: 

(H3) ^2 is the only sink in Pi and the unstable manifold of .^i is asymptotic to ^2- 

Therefore the union of the and their unstable manifolds forms a homoclinic cycle 
we label X. 

From the action of p~^, it is clear that P„ is spanned by a;„, Xi, . . . , xi. Using table[T] 
we deduce that Vi{c) is the and Vi{e) is the x^+i axis. The radial subspace of ^i 

has / dimensions, and so there are n — I — 2 transverse directions. Since has a one- 
dimensional unstable manifold, all the transverse and radial eigenvalues are negative. 

We label the contracting eigenvalue (that is, the eigenvalue in the direction spanned 
by Vj{c) at —c < 0, and the expanding eigenvalue (the eigenvalue in the direction 
spanned by Vj{e)), e > 0. We label the transverse eigenvalues — si, . . . , —Sn-1-2 where 
the eigenvector corresponding to the eigenvalue —Sk at is in the xi^j^k direction, 
taking subscripts modulo n as appropriate. 

Theorem 1 The homoclinic cycle X is asymptotically stable if 

n-l-2 
C + 5Z > 

If c + Yl]jZi~'^ < 6 'then the cycle is unstable. 

When the stability of the cycle changes at c + Y^^Zi "^ = 6, the cycle undergoes a 
resonance bifurcation. In this bifurcation a unique long-period periodic orbit bifurcates; 
the bifurcation may be supercritical (in which case the bifurcating periodic orbit is stable 
and exists for parameters for which X is unstable) or subcritical (and the bifurcating 
periodic orbit is unstable and exists for parameters for which X is stable). It turns 
out that the criticality of the resonance bifurcation is determined by the behaviour of 
trajectories in the 'global parts' of the flow near X (i.e. the parts of the flow away 
from neighbourhoods of the equilibria). In the second half of this paper we give an 
example of a system which contains a homoclinic cycle of the form described above, 
and we explicitly calculate these constants. For the example we give, we show that the 
bifurcation always occurs supercritically. The calculation also allows us to compute the 
period of the bifurcating stable periodic orbit. 

3. Proof of theorem [1] 

This section is devoted to the proof of theorem [TJ The proof is divided into three 
sections; first we construct return maps on Poincare sections in the standard way by 
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dividing the flow into 'local' and 'global' parts. The return maps describe the behaviour 
of trajectories asymptotically close to X. We then relate these maps to transition 
matrices, and finally deduce results on the eigenvalues of the transition matrices which 
enable us to prove the theorem. 

We remark that, although this result can be deduced from a theorem of Hofbauer 
and Sigmund [7] (see chapter 17, pages 220-232), we see substantial value in the proof 
presented here since it is more transparent for the situation at hand, and it makes the 
geometry of the invariant subspaces explicit. 



3.1. Poincare maps 

We define Poincare intersections about in the standard way: 

iff : {x„ = h'}, : {x^+i = h^}, 

for some small h > 0. We consider an initial point on if™: 

in,l / in,l in,l in,l in,l i 2\ 

^ ~ {'^1 ) • • • ) ; • • • 5 ^n-l5 J' 

where the Vj are radial coordinates which are zero at ^i. The 'time of flight' for the 
trajectory to reach is 

in,l\ 

Let the coordinates of the trajectory when it reaches H°'^*' be 

out,l _ out,l out,l out,l out,U 

then using the linearised flow at C,j, we have 

2r ■ / \ — 

out,l T i in.l / in,l 




OUL,_L 7 ^ 111.1 / 111,1 \ ■ -I J 

Vj =h - V- ' \x^^^ j , J = 1, . . . , / 

= h\ (10) 



OUtjl 

out,l . in,l / m,l 



out,i 7 i — : — ^ in,i / m,i 1 = ■ 7 I o 1 

' = h e Xj yx^^^ j ,j = / + 2,...,n-l 



This deflnes the local map near ^i. 

We now compute the global part of the Poincare map. The isotypic decomposition 
of with respect to the isotropy of W^[^\) is 

= © y,(si) © ... © \/i(s„_/_2) © Vl(c) 

where Vi(sj) = ((0, . . . , 0, x^+j+i, 0, . . . , 0)) is the eigenspace spanned by the eigenvalue 
Sj. Recall that dim Pi = Z + 1. The form of the global map from if^"* to a point x™'^ 
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on : {xi = /i^}, up to linear terms, is therefore: 



X 



in,2 



m,2 

,in,2 
1+2 



W2* 



Wl+1 





V / 



/ 

C21 



Q+11 




V 





C2l 





C21+1 








Q+iz Q+ii+i 

Cl+21+2 



\ 








/ out,l \ 

/ ^1 \ 



V. 



X 



out,l 

.out,l 
1+2 



■■0 ■■■ c 

where the Cij are 0(1) constants which depend on the global flow. By the invariance of 
the coordinate planes, Cjj > for I + 2 < j < n. In the section HI we find approximations 
for these cjj for a specific set of ODEs. The wj* are 0{h^) constants which define where 
the heteroclinic connection W'^{C,i) hits the plane H^. Note that x™'^ = h'^ by definition. 



and the radial coordinates on are 



in,2 



in,2 
^^l+l- 



We compose the local and global maps, and use the symmetry p to map points 



on onto H^, to give a return map 
with a point on H^, so 



on iJ™. That is, we associate the point p x 



l^m,2 



m,l^ 



(X 



m,2\ 



or 



/ ^1 \ 



X 



X 



m,l 

in,l 

l+l 
in,l 

1+2 



X 



in,2 



in,2 

'h+x 

Jn,2 
1+2 



/ in,2 
/ ^2 



in,2 

Jn,2 
«+2 



m,2 



The radial {vj) components are contracting, and do not affect the other components, 
so we only need consider an n — Z — 1 dimensional map ■?/' that describes the dynamics 
of the flow near the homoclinic cycle: 



in,l 
^1+2 



in,l 
^n-2 



/ ^/+2 \ 
m,2 



X 



Jn,2 
n— 1 



Vx-'V 



/ out,l \ 

/ l^/+22;i+2 \ 
out,l 



\ „out,l / 

\ ^n-^n / 



1+2 



+ 1 ) 

in.l 



Cn-lh 



°"-'-2 in,l /■ in,l 



V 



(12) 



where Cj = Cjj. For convenience, we relabel the coordinates and constants. We write 
p = n — I — 1, 

Vj = Xi+j, j = 1, . . . ,p, 
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to reach: 



= h'^Ci+j, for 1 < J < p - 1, Cp = h^^^--e^C„ 



Sj c 

aj = — , for 1 < j < p — 1, flp = -, 



(13) 



Vp-i Cp-iypVi"'' 



(14) 



The map ip has a fixed point at ?/i = . . . = ?/p = which corresponds to the 
homochnic cycle in the fiow. The map has a second fixed point, the yi coordinate of 
which satisfies 



where 



Note that 11%! C J = h^^^'^^ Uj=w^r ™^ ^^^^ P°i^* 

may correspond to a periodic 
orbit in the fiow, and we discuss the existence and stabihty of this fixed point further 
in section H] befow. 

In the next section we use transition matrix methods to determine the stabihty of 
the zero fixed point of ip and hence the stabihty of the homochnic cycle in the fiow. 



3.2. Transition matrices 

Let G be the set of mappings g : MP ^MP that have at lowest order the form 

g{xi, . . . , Xp) = [CiXi^^ X2^^ ■ ■ ■ Xp^^ , . . . , CpXi^ ■ ■ ■ Xp'''') 

for constants aij > and Ci non-zero. G is clearly closed under composition. We 
define the transition matrix [3], [19] of g to be the p x p real matrix M{g) with entries 
[M{g)]ij = aij. It is easily verified that if gi,g2 G G, then 

M{g2og,) = M{g2)M{gi). 

Any g & G has a fixed point at Xi = . . . = Xp = 0. The zero fixed point of the map g 
will be stable if all the row sums of M{g)^ diverge to +oo as — > oo. Conversely, if 
any of the row sums of M{g)'^ tends to 0, then the fixed point is unstable. 
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For the homochnic cycles described above, the transition matrix corresponding to 
the Poincare map f|T^ is: 



0-2 



1 






1 



ap_i 







■■■ 0\ 

... 

1 

0/ 



with the aj as in f|T3|) . We now recall the Perron-Frobenius theorem. We write M-v to be 
the product MijVj for a matrix M with components Mij and vector v with components 
Vj. The row sums of can be written as A^ ■ 1 where 1 is the vector (1,1,..., 1)-^. 

Definition 7 A real p x p matrix M is primitive if there exists an N > such that all 
entries of are strictly positive. 

By inspection, Ap is primitive, since aj > for all and so {ApY will have strictly 
positive entries. 

Theorem 2 (Perron-Frobenius) If M is apxp non-negative primitive matrix, then 
there exists a unique and simple positive eigenvalue Apf which is dominant in the sense 
that |A| < Apf for all other eigenvalues A of M. There exist right and left eigenvectors 
u, V > such that M ■ u = ApfU and • M = ApfV^. // u and v are normalised such 
that ■ u = 1 then 



hm (A/M^) 



T = uv' 



We now consider matrices of the form Ap, with Oj > 0, and in the remainder of this 
section, prove the following lemma. 

Lemma 1 IfYl^=i^j < ^> '^^s^ dominant eigenvalue Apf of Ap satisfies Apf <1. If 
X]j=i (^j > 1 then Apf > 1. 

Using the Perron-Frobenius theorem, we have 

lim X-f^A^ ■ 1 = T ■ 1 = u' > 0, 



where u' 



lim 



Af^oo 



If Apf < 1, then the row sums tend to zero, that is. 



Ap -1 = 0, and so the zero fixed point of the map ip is unstable. If Apf > 1, 
then Ap ■ 1 will thus be divergent to +oo, and the zero fixed point of the map will be 
stable. Since the stability of the fixed point in the map corresponds to the stability of 
the homoclinic cycle in the flow, substituting for the aj from ( |T3|) completes the proof 
of theorem [H 

We now prove lemma [H 



Claim 1 The eigenvalues A of Ap satisfy 

AP - (aiA^-i + aaA^-^ + . . . + a^.iA + Up) = 0. 



(15) 
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Proof. We will show by induction (on p) that 

Dp = det(Ap - XI) = (-If [A^' - (aiA^-^ + asA^"' + . . . + ap_iA + a^)] . 
Firstly, we have 

Di = ai-X = (-1)[A -ai]. 

Then 

ai-X 1 • • • 
02 -A 1 ••• 



(16) 



-XDk-i + i-lf-^ak. 



ak-i ■■■ -A 1 
ak ■■■ -A 
Assuming 0161) for p = A; — 1 we have 

Dk+i = -X{-lf [X'' - (aiA^-^ + asA'^-^ + . . . + ak-iX + ak)] + (-1) V+i 

= [A^+i - {a.X" + a2X'-' + . . . + a^A + a^+i)] . 

Hence f lTB]) is satisfied for all p > 1. □ 
The roots of equation (fTS!) are the eigenvalues of Ap. The Perron-Frobenius theorem 
states that there is at least one positive root, the next claim relates the size of the positive 
root(s) to the entries aj of Ap. 

Claim 2 Let X be a positive root of l[T5\) . If X]j=i < ^ then X < 1; if X]j=i > 1 
then A > 1. 

Proof. Consider the case X]j=i '^i < 1; suppose that A > 1: 

p 

AP = aiAP-i + aaAP-^ + . . . + a^.iA + < A^"^ ^ a^^ < A^"^ < A^. 

i=i 

This is a contradiction and hence A < 1. Similarly, for the case Yl^=i^j > 1' suppose 
that A < 1, and then 

p 

AP = aiAP-i + aaAP-^ + . . . + a^.iA + > A^"^ J] > AP"^ > Ap, 

i=i 

again a contradiction, and hence A > 1. □ 
Claim [2] is true for all positive roots A, therefore it is certainly true for the dominant 
eigenvalue Apf, completing the proof of lemma [H □ 



4. Bifurcation criticality computation 

In this section we give an example set of ODEs which contain a homoclinic cycle of the 
type discussed in section [2l For these equations, we are able to explicitly compute the 
shape of the heteroclinic connections between the equilibria, and use this to calculate the 
unknown constants Cj in the global part of the Poincare map f[T^ derived in section IXTl 
This allows us to compute the criticality of the resonance bifurcation and in addition, 
the period of the resulting bifurcating orbit. 
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4--1- System description 

The following equations satisfy conditions (HI) - (H3) and are equivariant under the 
action of V given in (I5l), ([6]) and (171). 

n 

ij — ~ 'y ^j -^i ~ ^-^j+l ~ Sn-3Xj_^2 — ■ ■ ■ — SiXj_2 + eXj_i), j = 1, . . . ,n (17) 

1=1 

for X = (xi, . . . , Xn) € M", c, e, Sj > 0, (j = 1, . . . , n — 3). In addition, in this example, 
1 = 1, and Q{x) is truncated at third order. We consider only n > A. The case n = 3 
is the example of Guckenheimer and Holmes [11] and it is well known that the system 
truncated at third order has a degenerate resonance bifurcation at c = e. 

The system f |T7j) has n equilibria on the coordinate axes we label the equilibrium 
has coordinates Xj = 1, Xk = 0, k 7^ j. Each equilibrium is hyperbolic, and the 
eigenvalues of the linearisation about are -2 (in the radial, that is, xj direction), 
— c (in the Xj^i direction), e (in Xj+i direction), and —Sk in the Xk+j+i direction 
{k = l,...,n — 3). Note that a homoclinic cycle between the equilibria C,j exists if 
—a and e are of opposite sign, but without loss of generality we consider only the case 
c, e > 0. Thus e is the expanding eigenvalue and c is the contracting eigenvalue. The 
—Sj are transverse eigenvalues; we assume —Sj < 0. By theorem [1], the homoclinic cycle 
is asymptotically stable if c + ^"=1 Sj > e. Otherwise it is unstable. 

In section [3TT] we constructed a Poincare map (fT2|) on a section H^, approximating 
the flow near the homoclinic cycle. This map has a non-trivial fixed point with X2- 
coordinate given by 

1 

(n \ i^s 

n (18) 
j=l+2 J 

where the Cj are defined in the map (fT2|) . 

This point corresponds to a periodic orbit in the flow, branching from the resonant 
homoclinic bifurcation, only if it is small as 5 ^ 1. Let C = ^J|"^;_^2 • ^'^^ C < 1 the 
fixed point is small if 5 < 1 (where the heteroclinic cycle is unstable), so the bifurcation 
is supercritical. Conversely for C > 1 the orbit exists in 6 > 1, and the bifurcation is 
subcritical. In the following, we write A = logC. Thus the bifurcation is supercritical 
if A < 0. 

4-2. The solution for the heteroclinic trajectory 

We now compute an approximation for the form of the heteroclinic trajectory near the 
resonance bifurcation, under the assumption that the difference between the expanding 
and contracting eigenvalues, e — c, is small. We write e = c + /3 and take |/3| <C 1. We 
consider the heteroclinic connection from ^1 to ^2 in the X1-X2 plane; all other connections 
are symmetry-related. The equations in the X1-X2 plane are: 

Xi = — [Xi + X2) — CX2), 

±2 = X2{1 - {xl + xl) + ex\). 
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For ease of computation, we make a change of variables, writing u = x1, v = x"^'- 

u = 2u{l-{u + v)-cv) (19) 

V = 2v{l- {u + v) + eu) (20) 

The heterochnic connection is a solution to ( fT9l) and (1201) with boundary conditions 
{u = 0,v = 1} and {v = 0,u = 1}. When c = e (/3 = 0), there exists an exact solution 
V = vo{u) = 1 — u. 

Recall that at resonance, 

c + s = e, 

where s = Yl^Zi ^j- Since e = c + at resonance we also have 
s = (3. 

Therefore, close to resonance, s is 0{(3), although it is important to note that it is not 
necessary that each sj is individually 0{(5). We will now look for approximate solutions 
to (HH]) and (I20D when /3 ^ 0. Write 

v = v^{u)+(3vi{u) + 0{(3''). (21) 

Expanding equations ( |T9l) and ( l20i) in powers of /?, we find at 0{j3): 



which gives 



v'liu) -vi{u) = --, 

CM 1 — M C 



In the following we also find it convenient to write the heteroclinic connection as a 
function u{v), that is, we write the connection in the form 

u{v) = l-v + (3ui{v) + 0{(3^). (23) 

By rearranging fl23|) and substituting into fl2T|) we find that ui{v) = vi{l — v) , that is 

^.c?. Constants in the global map 

We now use the approximated form of the homoclinic connection ( 1211) to compute the 
constants Cj in the map (fT2|) . From (fT2|) we have 

in,2 

Q = ^, j = 3,...,n. (25) 

Consider a trajectory close to the X1-X2 plane, so that all other coordinates are small. 
Approximations to the 0:3-,. . .,-Xn equations in (|T7|1 are therefore (ignoring higher order 
terms in the small coordinates): 

^3 = 3^3(1 — {^1 + ^2) ~ -^i^i + 6^2)5 (26) 
Xj = Xj{l - {xj + xl) - Sj_2x\ - Sj.^xl), j = 4, . . . , n - 1, (27) 

Xn = XniX — {x\ + x\) — Cx\ — Sn-zX^). (28) 
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We consider a trajectory which starts close to the heterochnic cycle. This trajectory 
will approximately follow the heteroclinic trajectory in the X1-X2 plane, so we integrate 
along the trajectory, from H°'^^ to H^. So from ( l26l) for instance, we find: 

/ — dx3 = / (1 - (xi + xj) - sixj + exj) dt. (29) 

J fjout X3 JHfut 

We can compute the left hand side of (!29|) : 

r -d^3 = log (xf'') - log {xT'') = log I ) = log((73), 

Jjjont X3 V / ' J 

(where the last equality follows from fl25l) ). and we find similar expressions for the other 
log(C'j) by integrating equations flTTI) and (1251) . 

As noted at the end of section 14.11 the criticality of the bifurcation depends on the 
sign of A, where 

(n \ n 

j=3 J j=3 

Summing the integrated forms of equations (126|) to (128|) we find, again writing xf = u, 

and X2 = V, 

A = y^\ogCj= {{n- 2)(1 - (u + v)) - cu-Y^ Sj{u + v) + ev\<lt 

U -^^r V U J 

where again the integral is taken along the heteroclinic connection in the X1-X2 plane. 
Note that this expression does not depend independently on the Sj] it only depends on 
the sum s = X]j=i ^j- 

We now rewrite the integral as an integral in u. Recall that the heteroclinic 
connection can be written as w(u) = 1 — u + Pvi{u) + 0(/5^), and e = c + /3. From (flQll . 
we have 

ii = -2cm(1 -u)- 2(3{c + l)uvi{u) + 0{(3'^), 



so 



u 2cu(l — u) 1(? u{\ — uy 
To determine the bounds on the integral, note that on H°'^^, x\ = h"^, so v = h, and 
the M-coordinate of the heteroclinic connection is u = 1 — h + Pui{h) + 0(/3^). On H^, 
u = h. 

Thus, expanding in powers of P (recall that s = 0(P)) we find: 

fHS" 1 
A= {{n — 2){1 — {u + v{u))) — cu — s{u + v{u)) + ev) —du 

Jfjout U 

-h 



[ {-cu + c{l-u)- s + f3 [-{n - 2 - c)vi{u) + (1 - u)]} 

Jl~h+0u^(h) 

du + 0{(3^ 



'l~h+f3ui{h) 

C + 1 Vi 



2c2 U(l 
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which, after rearranging terms and tidying up, becomes 

2A = r ^-i^- ^) d« _ ^ r 1 + s Z"'' — ^ du 

Ji-h+pm(h) u{l - u) Ji_f, cu Ji_h cu{l - u) 

^ (n — 2 — c)Vl(u) ,^ \C+1 Vx(u) , r^, n1\ 

1 1 P t I ^ s /-^ 1 1 ^ 
— du / — du -\ — / — I du 



I~h+I3ui{h) 



1 — U U C Ji_h C Ji^h \ — u 



log fi ^ 



+ logfl + ^^^j + — -^(log/i-log(l-/i)) 



Since the functions Ui and f i are (in principal) known through the integrals flM|) and fl22l) 
the above expression determines the sign of A, and hence yields the criticality of the 
resonance bifurcation. The remainder of this section is devoted to the computation of 
A which is perhaps surprisingly algebraically complicated. 

4-4- Expansion of integrals 

There are two types of terms in A which we cannot yet express explicitly: the terms 
involving ui{h) and the integrals involving Vi{u). We will compute both of these by 
calculating a power series expansion for Ui. The two integrals in A which we need to 
compute are: 

d„,./'"'^iiM^d». (31) 

Jh uyL-u) 

and 



u{l — u) w{l — w) 



'-'j!iM,,..r-Md... (32, 

h {i-uy Jh 



For convenience, we write q = 1/c, then 



w J Jo Vl-^. 

In the integrals (131!) and (1321) we only need to know Ui{w) for w < 1, so we can write 
the factors (1 — ^)~'' and (1 — as power series. In the following, (^) is a generalised 
binomial coefficient, that is, 

'q\_ r(g + l) 



^kj T{k + l)T{q - k + 1)' 
where T{q) is the usual Gamma function. Substituting the expansions into ui{w) gives: 



k=0 ^ ^ k=0 ^ 
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k( ^\ ^ „..q+k+l 



' ' k Jq+k+1 



oo 



qw Cnw"' (33) 



n=0 



where 



fe=0 



y g + A; + 1' 



5^(n-A;)( A; )g + A; + l 



fc=0 



Lemma 2 T/ie coefficient Cn in the expression for ui{w) given in ([22P, can be written 



as 



q 1 _ gr(n + l)r(g + l) 

n > 1 (35) 



g + r(g + n + 2) 

= (36) 
Proof. See appendix. 

We now evaluate the integrals (|3T1) and (l32l) . We first consider ( 13T|) . and writing 

'Ui(w) = X^fcLi Cfcw'^ we find: 



/ \ dw = q / Co— + > Cfcw'' 

oo 

= qco [\ogw]l''' + qJ2 [j 



°0 1 u 



k=l 
oo 

gco(-log/i) + g^| + 0(/i). 

fc=i 



In the appendix we show that 

oo 

Cfc q 



E 



and so, using this result together with ( l36l) . we have 

^dM; = —\ogh- +0{h). 37 

g + 1 (g + 1)^ 

Now we consider fl32|) . We expand both ui{w) and (1 — w)^^ to find: 

i-'-^-./-'(£.')(£..-)- 
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q / y^^dkw'^dw, 
•^^ k=0 



.fc=0 

oo 



k+1 



k=Q 



dk 



k+1 



+ 0{h), 



where we have defined = Yl'i=o^j- evaluate d^ in the appendix to find: 



E 

k=0 



dk 



k+1 q 



Hence 



l-h 



-dw = l + 0{h). 



'h w{l-w) 

Finally, we use our expansion of Ui to find Ui{h) = qh YlkLo ^kh^i and so 



-l3qco + Oiph) = -13-^ + Oiph) 



1 + q 



and 



log 



i-^-^)^om. 



{31 



(39) 



(40) 



Therefore, substituting ([37D, ([Ml), (I39D and (gOD into ([30]) we have 



2A = -f3—^ + q{2s - (3) \ogh - (3q{n - 1) + (3{q + 1) 
g + 1 

+0{f3h) + 0{(3^), 
= 2q{s-(3) \ogh + (3 
which simplifies dramatically to give 



q + 1 



log h) 



+ 



q + 1 



nq + q 



q+1 



+ 0{l3h) + Oils'" 



A = q{s- 13) \ogh 



[3nq 



+ 0{l3h) + 0{p"). 



By definition, as we approach the bifurcation point s — > /5, and so, evaluating A at the 
bifurcation point we obtain to leading order 

[3nq 



A 



< 0, 



(41) 



and hence the bifurcation is found to be supercritical for all values of n > 4 and q = 1/c. 
Note that although our calculation only computes A to leading order in /3, since the 
criticality only depends on the sign of A, the bifurcation will be supercritical whenever 
13 is sufficiently small. 
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5. Bifurcating periodic orbits 

In this section we discuss the evolution of the unique stable periodic orbit produced in 
the resonance bifurcation for the system considered in section |H 

First we use the calculation of the return map to derive an approximate expression 
for the period of the orbit. We emphasise that the analytic result contains no 
adjustable parameters. Referring back to the calculation of the return map, combining 
equations ([9]) for the time spent in a neighbourhood of an equilibrium on the cycle 
with (|T8|) which gives the leading-order relation between the location of the fixed point 
of the Poincare map and the coefficients Cj, we obtain the following expression for 
the period P of the periodic orbit, as usual neglecting the time spent moving between 
neighbourhoods of the equilibria on the cycle: 

P = ^T = -- log iC^/^^-'A, 
e ^ 

where 5 = (c + s)/e and C = YYj=i+2^j before. Substituting logC = A = —n(3q/2 
where /? = e — c > and q= 1/c this expression simplifies to yield 

p ^ n npq ^ r?{e-c) 

e1{\-S) 2c(e-c-s)' ^ ' 

This expression has been derived near the resonance bifurcation (i.e. |e — c — s| ^ 1), 
in the case that |e — c| ^ 1 and therefore is expected to be asymptotically correct in 
the limits of small s and small e — c. 

We compare the analytic result (l42l) with numerical integrations of (fT7|) . Our 
numerical integrations are carried out with a standard 4th order Runge-Kutta scheme, 
with the return times Pi,P2,... to the cross-section xi = 0.2 computed by linearly 
interpolating points on the trajectory that lie on either side of the cross-section. Since 
trajectories may be assumed to converge exponentially to the periodic orbit (discounting 
the possibility that it is nonhyperbolic) a highly accurate extrapolation of the true period 
of the orbit may be obtained by employing the ansatz Pk = P + 6oe~^^'^ and eliminating 
the coefficients (which are supposed to be constant) bo and bi by considering a set of 
three return times Pk, Pk+m and Pk+2m- This leads to the extrapolation formula 

p p p2 

p _ -T fc+2m-rfc fc+m ^^g-j 
Pk+2m ~ '^Pk+m + Pk 

which we find to give excellent results, even when the computed return times Pj used 
in (B3|) are far from the true (ultimate) period P. 

Figure [T] confirms the accuracy of the approximation fH2]) . showing the period P as 
a function of e — c — s for four different values of s in the case n = 4 (where there is only a 
single transverse eigenvalue). The excellent agreement for s = 0.01 (the solid curve and 
diamond symbols) and s = 0.1 (dashed curve and triangle symbols) shows that when 
< s <^ c the approximate expression fH2|) is extremely accurate. Figure [2] confirms 
the accuracy in the limit < s -C c when c is varied at fixed s; the data for c = 1 are 
most accurately approximated by (H2|) . Finally, figure [3] illustrates the dependence of P 
on n as being P - not a scaling that one might intuitively have proposed. For this 
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Figure 1. Orbit period P plotted against tlie eigenvalue combination e — c — s for 
the periodic orbit created in the resonance bifurcation at e — c — s = 0. Parameters 
are n — c — 0.5. Calculations are carried out by varying e for fixed values of s. 
Numerical data and the approximation (I42p are shown for four values of s, from bottom 
to top: s = 0.01 (0, solid line); s = 0.1 (A, dashed line); s — 0.2 (□, dash-dotted line); 
s — 0.5 (x, dash-triple-dotted line). The approximation is excellent for small s but 
systematically overestimates the period for larger s. 



figure we have for convenience taken s = 0.1 for n = 4, 5, 6 and achieved this by setting 
Sj = s/{n ~ 3) for all 1 < j < n — 3. We remark that the accuracy of (H2l) does not 
appear to vary significantly with either increasing or decreasing n. 

Further numerical investigations in the case that s > but allowing one, or more, 
of the Sj to be negative, so that the transverse eigenvalues, —Sj are positive, (i.e. when 
e — c — s < 0, the homoclinic cycle is only essentially asymptotically stable rather than 
asymptotically stable) show that a periodic orbit bifurcates ate — c — s = and has a 
period that is extremely close to that for the Sj = s/{n — 3) case as long as s remains 
much less than e. For example, in the case n = 5, s = si + S2 = 0.1, c = 0.5, e = 0.7 
we observed that P = 41.169 for Si = S2 = 0.05. On increasing si (and correspondingly 
decreasing S2 so that si + S2 = 0.1 always) we observed numerically that the period of 
the periodic orbit remained within P = 41.17±0.04 for 0.05 < si < 0.65 before reducing 
rapidly for larger si: P = 35.4 at si = 0.655 and P = 30.15 at si = 0.66. It seems 
reasonable to expect that the shape of the periodic orbit near the equilibria changes 
substantially when transverse eigenvalues become of the same order as the expanding 
eigenvalue e. We leave a detailed investigation of the dynamics to be the subject of 
future work. 
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Figure 2. Orbit period P plotted against the eigenvalue combination e — c — s for 
the periodic orbit created in the resonance bifurcation at e — c — s = 0. Calculations 
are carried out by varying e for fixed values of c. Parameters are n = 4, s = 0.1. 
Numerical data and the approximation ()42l) are shown for four values of c, from bottom 
to top: c — 2 (0, solid line); c = 1 (x, dash-triplc-dotted line); c = 0.5 (A, dashed 
line); c — 0.25 (□, dash-dotted line). The approximation is excellent for large c but 
systematically overestimates the period for smaller c. Note that the data for c = 0.5 
(A, dashed line) is shown with the same symbols and line style in figures [1] and [3l 



6. Discussion 

In this paper we have discussed resonance bifurcations in a class of vector fields 
containing robust homoclinic cycles that connect equilibria on a single group orbit, 
each lying in an Z-dimensional hyperplane in R" . Using the well-known return map 
technique, and the construction of transition matrices, we proved a general stability 
result for this class of homoclinic cycles in section [31 We then discussed the resonant 
bifurcation that occurs when c + s = e, where s is the sum of the relevant transverse 
eigenvalues. 

In the case that the equilibria lie on the axes in R'^ it is possible to make substantial 
further progress and compute the period of the bifurcating orbit by integrating along 
trajectories to explicitly compute the global parts of the return map. To the best of our 
knowledge, this is the first calculation of its kind. 

The calculation rests on the assumption that the sum s of the relevant transverse 
eigenvalues is small compared to the leading stable and unstable eigenvalues (c and e 
respectively). The details of the calculation become algebraically rather complicated, 
but only involve computing a number of integrals; this is achieved through power series 
expansions and summations. Unlike most return map calculations where the global 
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0.0001 0.0010 0.0100 0.1000 

e — c — s 



Figure 3. Orbit period P plotted against the eigenvalue combination e — c — s for 
the periodic orbit created in the resonance bifurcation at e — c — s = 0. Calculations 
are carried out by varying e for fixed values of s = 0.1 and c = 0.5, keeping all the Sj 
equal as discussed in the text. Numerical data and the approximation (j42|) are shown 
for three values of n, from bottom to top: rt = 4 (A, dashed line); n = 5 (0, solid 
line); n — Q (□, dash-dotted line). No trend in increasing or decreasing accuracy is 
observed as n varies. Note that the data for n = 4 (A, dashed line) is shown with the 
same symbols and line style in figures [T] and [21 

parts of the map introduce undetermined coefficients whose magnitudes are rarely 
known, our resuhs contain no free parameters and are seen to agree well with numerical 
computations of the period of the bifurcating orbit. 

It is of interest that the calculation shows that for this class of systems the resonance 
bifurcation is always supercritical. It would be of interest to investigate extending the 
present work to find a correspondingly simple class of examples where the criticality of 
the bifurcation depended non-trivially on parameters. For instance, in [20] we give an 
example of a cycle which has a resonant bifurcation, the criticality of which depends on 
parameters. It would also be interesting to see to what extent the analysis here could be 
applied to heteroclinic cycles in Lotka-Volterra systems, such as those studied, mainly 
numerically, by other authors [HI [15l [161 E! ■ 

Appendix 

In this appendix we provide the proof of lemma [21 and the computation of the sums 
Sfcli ^ X^fcLo Fk which are used in section in the computation of A. 
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Proof of lemmalE It is trivial to check that Cq = l/{q + 1). For n > 1, let 

n 

f{n) = J2Hn,k) 

where 

1 + g + A; - 1\ / q \fn + q + l 



q q + k + 1 \ K J \n — k J \ n 
that is, fin) is the right hand side of (134|) divided by the right hand side of ( 135!) . We 
will show that f{n) is independent of n, and specifically, that f{n) = 1, which proves 
the lemma. 

Note that F{n, k) = for k > n + 1. Further, let 

, . _ -k{l + k + q){k -n + q){k - l- n- nq) 
^ ~ qn{l + nf{l-k + n) 

for k 7^ n + and define G{n^ k) = R{n, k)F{n, k). Note that G{n,n + 1) is well defined, 
and G{n, k) = for k > n + 2. It can be verified that 

F{n + l,k)- Fin, k) = G{n, k + 1) - G{n, k) 

Summing over all k gives 

oo oo oo oo 

J2 F{n + l,k)-Y, F{n, k) =J2 ^i^, k + l)-J2 G{n, k) (A.l) 

k=0 k=0 k=0 k=0 

n+1 n n n.+l 

^F{n+l,k)-^F{n,k) =^G{n,k + I) -^G{n,k) - G{n,Q) (A.2) 

fc=0 fc=0 fc=0 k=l 

/(n + l)-/H = -G(n,0) = (A.3) 

so f{n) is independent of n. It is simple to check that /(I) = 1, completing the proof. 

□ 

The following lemma is used a number of times in what follows: 
Lemma 3 V ^'^^ + ^'^^^^ = ^ f 1 - imt^^] (A,4) 

Proof. We use induction on k. It is simple to check that (]A.4p holds for k = 1. Then 
assuming flA.4p is true for some A; > 1, we find for k + 1, 

^^ r(a; + i)r(j + i) Ji^ r(x + i)r(j + 1) r(a; + i)r(fc + 2) 

^ r(a; + j + l) ^ r(x + j + l) r(x + A; + 2) 

r(x)r(A; + 2)\ xT{x)T{k + 2) 



X 


X 




1 




X 




X 




1 




X 




X 




1 



1 



T{x + k + l)J (x + fc + l)r(x + A; + l) 
r(x)r(A; + 2) / k + 2 



T{x + k + l) \x + k + 1 
T{x)T{k + 3) 
r{x + k + 2) 



□ 
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Corollary 1 For x > 1, 

r(a; + l)r(j + 1) X 



T{x + j + l) x-1 
Proof. The required result is equivalent to showing that 

fc^oo T{x + k + 2) 

To show this we write the combination of Gamma functions as a Beta function: 

T{x)T{k + 3) _ ^ + 2, x) = (x - l)B{k + 3, X - 1), (A.5) 



r(x + A; + 2) 
where we define the usual Beta function 

pi poo q-l 

B(p,q)= / tP-\l-ty-Ut= / r^ds, (A.6) 

^ ^ Jo Jo {i + sy+1 ^ ^ 

making the change of variable 1 + s = t~^, see [21]. The last equality in flA.51) can be 
derived in a straightforward manner by integrating by parts the first integral expression 
in flA.6p for B{p, q). Then we have 

(x-l)B{k + 3,x-l) = ^^—-L^^ds< {x-l){l + s)-'-Us, 

{x-i){i + sy'-'^ 



< - 



k + 3 



s=0 



X — I 

A; + 3' 



which clearly tends to zero as /c ^ cxd for any fixed x > 1. □ 
We are now able to compute the following two sums which are used in section 14.41 
Firstly, consider 

^c, ^_^gr(fc + i)r(g + i) 

k=l 



k kV{q + k + 2) ' 



q y. T{k)T{q + 3) _ q r(A; + l)r(g + 3) 



(g + 2)(g + l)^r(g + /c + 2) (g + 2)(g+l)^ r(g + A; + 3) 
q q + 2 q 



(g + 2)(g + l)g + l (q + l)^' 

Now consider 

k 
j=0 

Ar(j + i)r(g + i) 

= Co-q2^ 



Tiq + j + 2) ' 
1 , gr(g+l) Ji^r(j + l)r(g + l) 



1 + g r(g + 2) r(g + j + 2) 

J_,^ 9 r(j + l)r(g + 2) 

1 + q 1 + q g + r(g + j + 2) ' 
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^ q + _ r(g + l)r(A; + 2) 



(g + 1) g V r(g + A; + 2) 
^ r(g + l)r(A; + 2) 
r(g + A; + 2) ' 
so that the second sum can be computed to be 

^k+l ^ Tiq + k + 2) ' 

fc=0 A;=0 ^ 



oo 



r(g + 2)r(A: + 1^ 

5 + 1^^ r(g + /c + 2) 

g + 1 g g 
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